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I.  Introduction 

The  e-f-fect  o-f  suspended  particles  in  a  water  medium  on  the 
transport  o-f  light  is  normally  dichotomized  into  two  processes. 
Namely,  the  dual  e-f-fects  o-f  absorption  and  scattering. 
Absorption  simply  acts  to  remove  photons  which  might  have 
contributed  to  image  -formation.  Such  attrition  e-f-fects  might  be 
overcome  -for  a  given  system  by  increasing  the  souce  -flux  or  by 
using  a  more  sensitive  receiver.  Scattering,  however,  produces 
two  more  pernicious  e-f-fects.  Scattering  can  produce  a  -foreground 
veiling  glare  which  reduces  the  e-f-fective  contrast  o-f  the  target 
with  respect  to  the  background,  and  by  causing  a  redistribution 
o-f  the  trajectories  o-f  image  -forming  photons,  scrambles  the 
in-formation  content  originally  presented  by  the  target  and 
background. 

The  modulation  trans-fer  -function  (MTF)  is  used  to 


characterize  the  e-ffects  the  medium  has  on  the  passage  of  image 
forming  light.  It  is  not  surprizing  that  values  of  optical 
parameters  which  describe  the  scattering  effects  of  the  water 
medium  are  significant  variables  in  the  MTF  expression.  Wells^-^ 
has  derived  a  transformation  which  converts  the  medium  volume 
scattering  function,  P<6)  ,  into  a  medium  MTF  expression  via  the 
use  of  a  'decay'  function.  The  MTF  expression  so  derived 
requires  the  value  of  the  total  (volume)  scattering  coefficient, 
b.  The  value  of  b  i s  related  to  the  volume  scattering  function 
by  the  functionality 

2>>2ff  r”p  (0)  sinBdD  ( 1  > 

^0 

stated  here  in  spherical  coordinates  and  assuming  azimuthal 
symmetry.  Normally,  b  is  not  measured  directly  due  to  either 
lack  of  appropriate  instrumentation  or  simply  the  time  consuming 
nature  of  the  task.  Rather,  the  value  of  b  is  computed  by  using 
other  more  easily  determined  optical  coefficients. 
Prei sendorf er^ ,  Honey*,  Wilson”,  Timofeyeva*,  and  Morel ^  provide 
some  specific  relationships  which  can  be  used  to  determine  b  for 
a  given  water  mass. 

The  objective  of  this  paper  is  to  utilize  Wells'  results  to 
compute  the  medium  decay  function,  D(D ,  (from  which  Is  easily 
determined  the  germane  MTF  for  a  stated  range,  R>  using  both  the 
value  of  b  computed  with  the  above  mentioned  relationships  and 
experimentally  determined  values  of  b.  The  resulting  decay 
function  expressions  will  be  compared  to  provide  a  measure  of  the 


errors  introduced  by  using  the  various  estimation  methods. 
Finally,  comments  will  be  made  on  which  model  appears  to  provide 
the  most  accurate  estimation  o-f  b. 

I I . Background  and  Methods 

The  medium  MTF  is  the  sine  wave  -frequency  amplitude  response 
of  the  water  path.  Wells'  transf ormati on ^ =  of  the  volume 
scattering  function  into  a  medium  MTF  assumes  the  small  angle 
approximation  and  that  the  suspended  particles  are  larger  than 
the  optical  wavelength. 

Wells'  expression  for  the  medium  MTF  can  be  written  in  the 

f  orm= 

T(r,J?)-e3q>I-I?(r)J?)  <2) 

where  R  is  the  range  of  Interest  and  D(I^  is  the  'decay'  function 
with  r  the  spatial  frequency  in  cycles/radian.  The  decay 
function  may  be  written  as^ 

D(r)-i>-o(r)  <3) 

where  Q(D  provides  for  a  recovery  of  power  at  lower  spatial 
frequencies  where  the  scattering  is  too  small  to  create  blur  for 
target  constituents  whose  relative  dimensions  are  large^.  Here 
it  is  assumed  that  the  scattering  function  is  highly  peaked  in 
the  forward  direction.  It  is  important  to  notice  that  the  decay 
function  is  Independent  of  the  range,  R.  Q(D  is  defined  by*- 


where 


i?<r)-2xj^*u(0)  Jo(«)d9 


(4) 


u(d)-f*p(8)ci5  (5) 

Jo 

and  Jai(a)  is  a  zeroth  order  Bessel  function  with  «  =  2sr0. 

It  can  be  shown that  QCF)  monotonical  ly  decreases  as 
spatial  frequency  increases  (i.e.  as  the  potential  for  image 
detail  increases),  and,  in  the  limit  of  high  spatial  frequency, 
Q(D  approaches  zero.  In  this  region,  D(D  approaches  b.  Hence, 
the  diminuition  of  fineness  of  resolution  of  the  image  forming 
light  over  a  range,  R,  is  due  to  an  exponenti al -type  decay  of  the 
relevant  spatial  frequencies  which  is  controlled  by  the  value  of 
b. 

Table  1  provides  the  relationships  used  to  estimate  b 
knowing  other  optical  parameter  values  for  the  water  mass  in 
question. 

Table  II. 


Name 

Rel ati onshi p 

Parameters  Used 

Honey 

b=6/5*(c-k) 

c,k 

Wi 1  son 

b*l/0.85*(c-k) 

c,k 

Prei sendorf er 

b=c-3/4*k 

c,k 

Timofeyeva 

k/c=(0.  19*(l-b/c) 

c,k 

Morel 

b=0.30*<Chl  >«••«»= 

Chi 

In  Table  II,  c  represents  the  beam  attenuation  coe-f -f  i  ci  ent ,  k  the 
dif-fuse  attenuation  coefficient,  and  Chi  the  concentration  of 
Chlorophyll  a. 

The  scattering  data  used  in  the  computations  were  collected 
in  a  large  research  tank  facility  during  laser  backscatter 
experiments  conducted  at  the  Naval  Oceanographic  and  Atmospheric 
Research  Laboratory®  (NOARL) .  The  tank  water  was  cleaned  by 
continuous  pumping,  filtering,  and  sedimentation.  Then, 
beginning  from  the  initial  clear  water  state,  the  water  was  made 
increasingly  more  turbid  in  a  controlled  sequence’’.  Ancillary 
supporting  optical  measurements  were  made  among  which  were  water 
transmission,  diffuse  attenuation  coefficient,  and  volume 
scattering  function.  The  volume  scattering  function  measurements 
were  fit  by  a  modified  Heyney-Greenstei n  function*"*  so  that 
values  of  the  total  scattering  coefficient  could  be  computed  for 
the  different  water  mass  turbidities'’.  The  relevant  wavelength 
for  the  tank  data  used  in  this  work  is  528nm. 

III.  Results  and  Discussion 

As  shown  in  equation  (2),  the  MTF  can  be  determined  by  the 
negative  exponential  of  the  product  of  the  decay  function  and  the 
range.  Hence,  to  make  relevant  calculations  independent  of  the 
range,  the  decay  function  will  be  computed  for  the  germane  b 
values.  The  MTF  can  then  be  formed  easily  by  use  of  equation 


(2) 


Figure  1.  Wells'  decay  function  as  a  function  of  spatial 
frequency  and  value  of  b  used. 


Figures  1-3  provide  the  decay  function  for  different 
experimental  and  estimated  values  of  b.  Here,  it  is  understood 
that  the  values  of  the  decay  function  have  not  been  divided  by 
the  arguement  value  of  the  Bessel  function.  Figure  1  provides  5 
curves.  The  solid  line  is  the  decay  function  computed  with  the 
experimentally  determined  b.  The  other  curves  shown  provide  the 
decay  curves  for  b  estimated  using  the  methods  given  in  the 
legend;  and  similarly  for  Figures  2  and  3. 


Figure  2.  Wells'decay  -function  for  an  intermediate  value  of 
b  from  test  tank  data. 


The  average  error  is  given  in  Table  2. 

Table  2 . 


Method 

Error 

T i mof eyeva 

4 1 . 77. 

Honey 

7.47. 

Prei sendorf er 

8.67. 

Wi 1  son 

6.97. 

Morel 

40.  17. 

The  average  error  was  computed  by  taking  the  mean  o-f  the 
di-f Terence  between  and  estimators  values  oT  the  decay  -function 
and  the  experimentally  determined  b  decay  -function. 

Essentially,  the  curves  indicate  how  accurate  the  methods 
•for  estimating  b  are  when  compared  to  actual  b  values.  Figure  4 
provides  a  plot  of  b  vs  c  values  from  the  research  tank  data. 
Here,  it  is  seen  that  the  Honey/Wilson  (since  these  two  are  very 


Figure  3.  Mells'  decay  function  for  a  high  value  of  b. 
provides  the  most  accurate  estimate  of  b. 

However,  the  values  of  b  for  the  tank  were  relatively 
high.  Perhaps  in  real  world  waters,  with  lower  values  of 
scattering,  a  different  method  of  estimating  the  value  of  b  would 
be  more  successful.  To  check  this,  a  NOARL  dataset  for  the 


Sargasso  Sea  region  was  used  for  which  the  basic  quartet  of 


irradiances  were  measured.  Figure  5  exhibits  the  results  o-f  the 


calculation.  It  is  seen  in  these  very  clear,  oligotrophic 
waters,  over  a  small  range  o-f  c  and  b  values  determined  by  the 
dataset,  the  Prei sendor-f er  relationship  appears  to  provide  the 
most  accurate  estimate  o-f  b  overall.  This  is  closely  -followed  by 
the  Wilson  relation. 
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H 
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-•-to  DAIA  -•.TIMOFEYEVA  HONEY 
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Figure  4.  Research  tank  data  results  -for  b  values. 

In  our  computations,  it  is  assumed  that  the  water  path  o-f 
interest  is  horizontal  (that  is,  a  plane  parallel  strati-fied 
medium  is  assumed).  l-f  slant  or  vertical  paths  are  considered, 
then  each  layer  o-f  the  medium  represents  a  di-f-ferent  set  of 
optical  conditions  —  much  like  a  new  lens  element  of  an  optical 
system.  Thus,  for  these  cases,  an  MTF  may  be  computed  for  each 
layer  and  the  result  cascaded  to  provide  an  overall  medium  MTF. 
The  b  values  to  be  used  for  such  medium  MTF's  must  then  be 


pro-files  the  b  value.  That  is,  values  o-f  b  as  a  function  of 
depth.  Using  the  Honey/Wilson  or  Prei sendorf er  approach, 
profiles  of  k  and  c  are  needed.  These  may  be  obtained  with 
standard  optical  instrumentation.  However,  the  Morel  relation 
requires  a  knowledge  of  the  Chlorophyll  a  profile.  Often,  this 
is  obtained  with  a  fluorometer  which  is  calibrated  to  read  the 
Chlorophyll  a  values  in  the  water.  Kitchen  and  Zaneveld^^  have 
recently  shown  that  the  Chlorophyll  maxima  do  nc  .  necessarily 
coincide  with  all  scattering  present  due  to  the  possibility  of 
seston  layer  scattering  at  other  depths.  Thus,  if  the  Morel 
relationship  is  used,  some  means  of  determining  additional 
scattering  layers  independent  of  Chlorophyll  concentration 
appears  needed. 

IV.  Summary 

Decay  functions  were  calculated  for  both  experimentally 
determined  and  estimated  values  of  b.  The  relevant  MTF  may  be 
calculated  directly  from  the  decay  function  curve  given  the 
range,  R,  of  interest  (see  equation  (2)).  A  table  listing 
average  error  over  the  b  values  used  was  given. 

Of  the  relationships  tried  (Table  I),  the  Honey/Wilson 
relation  proved  to  give  the  most  accurate  estimation  of  b  for  the 
research  tank  data.  However,  for  the  Sargasso  Sea  dataset,  the 


relation  of  Prei sendorf er  gave  the  better  rendition  of  b. 
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Figure  5.  Values  O'f  b  -from  a  Sargasso  Sea  dataset. 
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